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Abstract 

We introduce an alternative discretization for static quarks on the lattice retaining the 
0(a)-improvement properties of the Eichten-Hill action. In this formulation, statis- 
tical fluctuations are reduced by a factor which grows exponentially with Euclidean 
time, xq. For the first time, B-meson correlation functions are computed with good 
statistical precision in the static approximation for xo > 1 fm. At lattice spacings 
a ~ 0.1 fm, 0.08 fm, 0.07 fm, the B s -meson decay constant is determined in the com- 
bined static and quenched approximation. A correction due to the finite mass of the 
b-quark is estimated by interpolating between the static result and a recent determina- 
tion of Fy) s . 
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1. B-physics matrix elements such as the B-meson decay constant Fb are obtained 
from lattice correlation functions at large Euclidean time. Considerable interest 
lies in the treatment of the b-quark in the leading order of HQET, the static 
approximation [1,2]: in this framework non-perturbative renormalization can be 
performed, the continuum limit exists and also corrections can in principle 

be taken into account [3-6]. 

Progress along this line has been hampered by large statistical errors in the 
static approximation. In particular it has been observed [7] that the errors of a 
B-meson correlation function roughly grow as 

noise 

-Rns = 7 oc exp (x Q A) , A = £ stat - m n , (1) 

signal 

where -E s t a t is the ground state energy of a B-meson in the static approximation 
with the Eichten-Hill action 1 , 

S h EH = a A Y^^)D^{x), (2) 

X 

DoMx) = -f/ t (x-a6,0)Vh(x-a6)], (3) 

a 

for the static quark [2]. Eq. ((TJ) is problematic because the requirement -Rns "C 1 
is satisfied only for xq of the order of A -1 and this time interval shrinks rapidly to 
zero in the continuum limit a —>■ where -E s t a t ~ &i x <7o/° with some number e\. In 
the attempt to eliminate the discretization errors by reducing the lattice spacing, 
a, one is then limited more and more by unwanted contaminations by higher 
energy states and it has been very difficult to compute matrix elements in the static 
approximation [1,8-11]. Since the exponent in eq. ((TJ) is dominated by a divergent 
term, it is plausible that one may reduce it by changing the discretization. Here 
we will demonstrate that this is indeed possible while remaining with roughly the 
same discretization errors. 

In [3] it has been shown that energy differences computed with the action 
eq. © are 0(a)-improved if the relativistic sector (light quarks and gluons) is 
0(a)-improved. Furthermore, apart from the usual mass dependent factor, 1 + 
6^ at am q , the static axial current, 

Af t (a;)=V; i (x)7o75^h(x), (4) 
is on-shell 0(a)-improved after adding only one correction term, 

(Af* ) = + adpSAg* , 5A^(x) = Ux) ^ 3 + ^ W ) • (5) 

^or a more precise definition of the theory and for any unexplained notation we refer to [3]. 
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We want to retain these properties of the theory. They are guaranteed if the 
lattice Lagrangian is invariant under the following symmetry transformations (we 
do not list the usual ones such as parity and cubic invariance) [3]. 



i) Heavy quark spin symmetry: 

ip h — > Vip h , Vh — > ^bV' 1 > witn V = exp(-i(piei jk cr jk ) 



ii) Local conservation of heavy quark flavor number: 



-ir?(x) 



(6) 



(7) 



Keeping these symmetries intact, there is little freedom to modify the action. 
We may, however, alter the way the gauge fields enter the discretized covariant 
derivative, Do. To this end we choose 



D ip h (x) = - [ip h (x) - W\x - ab, 0) VhO - a0)\ 

CL 



(8) 



with W(x, 0) a generalized gauge parallel transporter with the gauge transforma- 
tion properties of U (x, 0). In particular we take W(x, 0) to be a function of the link 
variables in the neighborhood of x, which is invariant under spatial cubic rotations 
and does have the correct classical continuum limit such that D Q = do+A Q + 0(a 2 ). 
This is enough to ensure that the universality class as well as 0(a)-improvement 
are unchanged in comparison to eq. (JHJ)- Since we expect the size of remaining 
higher order lattice artifacts to be moderate if one keeps the action rather local, 
we here consider only choices where W(x, 0) is a function of gauge fields in the 
immediate neighborhood of x, x + a0. We choose 



W s {x,0) 

Wa(x,0) 
W nYP (x,0) 



V(x,0) 

V(x,0), 
Vkyp(x,0) 



|+(itryt( X)0 )F(x,0) 



1/2 



(9) 
(10) 
(11) 



where 



1 3 

^,0) = -^ [ U(x,j)U(x + a],0)U^x + a6,j) 



3=1 



+ rf(x-aj,j)U{x-aj,Q)U{x + a6-aj,j) ], (12) 



and where the so-called HYP-link, Vuyp(x, 0), is a function of the gauge links 
located within a hypercube [12, 13]. In the latter case we take the parameters 
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cti = 0.75, «2 = 0.6, 0:3 = 0.3 [12]. The choices (JOJ) - ((Hj) will be motivated further 
in [14]. It is worth pointing out that a covariant derivative of the general type 
used above has first been introduced in [15]. In this reference it was considered for 
the Kogut-Susskind action for relativistic quarks and with a different motivation. 

2. Next we have to study the scaling behaviour of observables computed with 
the actions S^, S^, S^ YF which are obtained by inserting W$, Ha, W^hyp into 
eqs. (JBI) and (J2J). In [14] this scaling behaviour is analyzed in depth for various 
observables and various choices for the static action in perturbation theory and 
non-perturbatively. Here we will present only one example. The necessity of such 
an investigation can be underlined by the following consideration. 

The static potential can be seen as an energy for a static quark with action 
5h and an antiquark with the corresponding [16]. Hence, the static force is one 
indicator for the scaling behavior of these actions. In [13], rather large a 2 -effects 
have been seen in the short-distance force for S^ YP and S^ YP . 

One may therefore worry about large a-effects, in particular in correlation 
functions of the static-light axial current, where static and light quarks propagate 
also close to each other. With the new actions, A S Q at is 0(a) -improved once [14] 

-0.08237£ 2 + OG7 4 ), forS h = S h EH , (13) 
0.0072(4) g 2 + O(g*), for S h = Sls£, (14) 
0.0385(37) g 2 + O(g 4 ), for S h = S h HYP , (15) 

is set in eq. The improvement coefficient 6^ at is set to its tree-level value 
^stat _ in this work. 

We consider now a step scaling function, E^ at , which gives the change of the 
renormalized static axial current in a Schrodinger functional (SF) scheme [17], 
when the renormalization scale is changed from fi = 1/L to fi = 1/(2L). Its 
continuum limit is known for a few values of L [17]. This quantity is thus a good 
observable to search for a-effects. In Fig. we show E^ at (3.48, a/L), where the 
first argument parameterizes L in terms of the SF-coupling g 2 (L) = 3.48. 0(a)- 
improvement is employed as in [17] but we consider the different actions for the 
static quark introduced above. All of them lead to S^ at (3. 48, a/L) at finite a/L 
differing from the continuum limit by about the same amount. Supported also by 
further such studies [14], we conclude that within the set of actions studied none 
is particularly distinguished by its scaling behavior. 

3. Let us now demonstrate that the statistical errors at large Euclidean time 
are reduced by the choices eqs. © - (fTTj). As a B-meson correlation function we 



stat 
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stat 
C A 

stat 
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Figure 1: The step scaling function 
Results for S® B were extrapolated to S^ at (3.48, 0) [17] (•). In all cases c^ at from 1-loop pertur- 
bation theory is used, which is sufficient since £^ at (3.48, a /L) does not depend very sensitively 
on this improvement coefficient. For S£, and S^ YP points have been displaced on the 
horizontal axis for clarity. 



choose 

/r (*o, ") = -U(ArUx) 0{u)) , 0{u) = ^ J2 Ch(y)7 5 ^(y - z)Ci(z) , 

(16) 

defined in the Schrodinger functional with T = 3L/2, L/a = 24, (3 = 6/g% = 6.2 
and a vanishing background field [3]. Here, as a novelty compared to previous 
applications, a wave function u>(x) is introduced to construct an interpolating B- 
meson field in terms of the boundary quark fields (j and ( h . It may be exploited 
to reduce the contribution of excited B-meson states to the correlation function, 
but this does not concern us yet. At the moment we simply consider cj(x) = 1 
and form the ratio i?NS, eq. (JIJ, for the different actions. From now on we set the 
light quark mass to the strange quark mass, taken from [18] following exactly [19] 
concerning the technical details. 2 Figure 121 shows that in all cases R^s grows 
exponentially with xq. For the Eichten-Hill action, also the effective coefficient 
A, describing the growth for xq = lfm — 2fm, is roughly given by E stSut — 
in agreement with eq. while for the other actions this is not the case. Most 
importantly for the other actions, A is reduced by a factor around 4, and with 
the statistics in our example a distance of xo ~ 2 fm is reached with S^ YP if one 
requires i?NS < 2%. The actions S^, Sf behave only slightly worse. 

2 Of course, these details matter only before taking the continuum limit. 
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Figure 2: The ratio Rns, eq. Q, for the correlation function /j[ tat for a statistics of 2500 
measurements. Circles refer to Sjp 1 while squares and triangles to Sf and S^ YF , respectively, 
behaves like Sf . Physical units are set by using ro = 0.5 fm [20,21]. 



4. This reduction of statistical errors enables us to choose u;(x) such that a long 
and precise plateau is visible in the effective energy, 



E eS (x , u) = In [/f*(a; - a, u)/ff & \x Q + a,u)] /(2a) 



(17) 



as shown in Fig. |31 Neither position nor length of the plateau depend sensitively 
on the details of a;, as long as it is chosen such that the first excited state in the 
B-meson channel is canceled to a good approximation. For the figure as well as 
for the following, we have chosen u G {Qi, f^} with 



57i(r) 
w 4 (x) 



uji + ato 3 , Vt 2 = UJ2 + ol'uja , 
N- 1 ^SJiflx-nLl), £ = 1, 2, 3 



ngZ 3 

r -3/2 e - r /a ^ = ^3/2 e -r/(2oo) ; = r ~5/2 f e -r/(2oo) 



(19) 



where ao = 0.1863r and the (dimensionless) coefficients Ni are chosen such that 
a 3 ^ x tu 2 (x) = 1. The B-meson decay constant is then obtained from the renor- 
malization group invariant matrix element [22] 



$ rgi(zo) = ~Z RG1 (1 + fc^am,) 2L 3 / 2 
of the static axial current, where 



fstat/™ \ 

JA \ X 0) Jx -T'/2)E eS (xo) 



(20) 



h{T,u) = -i(C( w )0( w )>, = f I ^Ci'(y)7 5 ^(y-z)C h '(z). (21) 
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10 20 30 x /a 

Figure 3: Effective energies for wave functions Qi (open symbols) and SI2 (filled symbols) using 
S£ YP (circles) and Sf (triangles). Results refer to a 24 3 x 36 lattice, f3 = 6.2. 

The renormalization factor, Zrgi, relates the bare matrix element to the renor- 
malization group invariant one [17]. Its regularization dependent part is computed 
exactly as in that reference, but for the new actions. In Tabled we give results for 
$rgi(^o) for three values of the lattice spacing and selected choices of T,T',Xo, 
highlighting what we selected for further analysis. These numbers do not change 
significantly if we vary the improvement coefficients c^ at and 6^ at , which are known 
only in perturbation theory, by factors of two. We thus extrapolate our results 
quadratically in the lattice spacing and arrive at our estimate for the continuum 
limit 

r 3/2 $RGi = 1-78(13) . (22) 



5. The result eq. (|22|) may be used to compute Fb s by taking account of the 
mass dependent function [17] Cps(Mb/A^jg) = -Fba/^b/^rgi = 1-22(3), evaluated 
using the 3-loop anomalous dimension [23] and the associated matching coefficient 
between HQET and QCD [24]. Mb denotes the renormalization group invariant b- 
quark mass [4,5]. With this we arrive at tqF^ = 0.57(4). A correction due to the 
finite mass of the b-quark can be computed by connecting the static result eq. (|22|) 
and 

r ° 2 c F ui7l' ) = l33<7) < 23) 
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16 


24 


20 


10 


1.851(17) 
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-0.200 


6.1 


0.079 


24 


30 


30 


15 
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0.022 
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30 


12 
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0.756 


0.022 


6.2 
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24 


36 


36 


18 


1.724(78) 


1.760(75) 


0.351 


-0.176 


6.2 


0.068 


24 


36 


36 


15 


1.726(78) 


1.763(76) 


0.351 


-0.176 



Table 1: Decay constant in static approximation. 



by a linear interpolation in the inverse meson mass. Here we have used recent com- 
putations of the D s -meson decay constant [25] and of the charm quark mass [19]. 
In this way we obtain 

r F Bs = 0.52(3) -> F Bs = 205(12) MeV with r = 0.5 fm . (24) 

Conservatively, we may attribute an additional ~ 5 % uncertainty to the fit ansatz 
used, but our personal estimate is that this error is significantly smaller and it 
will soon be quantified [26]. 

One should remember that eq. ()24j) refers to the quenched approximation and 
as in [25] a 12 % scale ambiguity may be estimated from the slope of the linear 
interpolation. 

6. An interesting point is that the potential in full QCD may be computed 
replacing the time-like links in the Wilson loop (or Polyakov loops) by the different 
Wi introduced above. In particular the "HYP-link potential" [13] may be used. 
Depending on which Wi is chosen, the static potentials differ from each other, 
but all of them approach the continuum limit with 0(a 2 ) corrections if the action 
used for the dynamical fermions is 0(a)-improved. This property follows from 
the considerations of [16] applied to the static actions introduced above, which 
satisfy all the necessary requirements. This virtue of e.g. the HYP-link potential 
was not obvious before. Using it, better precision can be reached and some signs 
of string breaking [27] may become visible. 

7. To summarize, we have shown that a modification of the Eichten-Hill static ac- 
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tion can be found which keeps lattice artifacts in heavy-light correlation functions 
moderate but reduces statistical errors to a level making the region xq > 1.5 fm 
accessible. Furthermore, the new action can be used without change for dynam- 
ical fermions and also to compute the static potential with dynamical fermions. 
As a demonstration of the usefulness of this reduction of statistical errors, we have 
computed _Fb s in the quenched approximation, by joining the continuum limit of 
the static approximation estimated with the new action with the previously deter- 
mined continuum limit of Fd b by means of a linear interpolation. This procedure 
can systematically be improved by computing 1) the mass dependence around m c , 
2) the 1/m corrections to the static limit and 3) repeating the whole analysis with 
dynamical fermions. Work along these lines is in progress and a more detailed 
investigation of the properties of various static quark actions is in preparation. 
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